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ABSTRACT 

We present an analysis of the X-ray properties of 35 submm galaxies (SMGs) in the 
Chandra Deep Field North region. Using a sample of robust 850 /im-selected galaxies, 
with sub-arcsecond positions from Spitzer and/or radio counterparts, we find 16 ob- 
jects (45 ±8 per cent) with significant X-ray detections in the 2-Ms Chandra data. 6 of 
these SMGs (~ 17 ±6 per cent) have measured X-ray luminosities or upper limits con- 
sistent with those expected based on the far-infrared or radio-derived star formation 
rate, and hence with the X-rays coming solely from star-forming processes. Extrap- 
olating observed X-ray/star formation rate relations to the luminosity of the SMGs 
we find that the X-ray derived star formation rates are typically in the range 200- 
2000 Moyr- 1 . In another 7 sources (20 ± 7 per cent of the SMGs) a dominant AGN 
contribution to the X-ray emission is required, while in 3 more it is unclear whether 
stellar process or accretion are responsible for the X-rays. Stacking of the X-ray un- 
detected SMGs reveals a highly significant (7a) detection. Under the assumption that 
the stacked X-ray is due to star formation this gives an average X-ray derived star 
formation rate of around 150MQyr _1 . We deduce that the AGN fraction in SMGs 
based on X-ray observations is 20 — 29 ± 7 per cent, which is towards the lower limit 
of previous estimates. Spectral analysis shows that in general the SMGs are not heav- 
ily obscured in the X-ray but most of the SMGs classfied as AGN show absorption 
with column densities in excess of 10 22 cnr 2 . Of the secure AGN, the bolometric lu- 
minosity appears to be dominated by the AGN in only 3 cases. In around 85 per 
cent of the SMGs, the X-ray spectrum effectively rules out an AGN contribution that 
dominates the bolometric emission, even if the AGN is Compton thick. The evidence 
therefore suggests that intense star formation accounts for both the far-infrared and 
X-ray emission in most SMGs. We argue that, rather than having an especially high 
AGN fraction or duty cycle, SMGs have a high X-ray detection rate at very faint 
fluxes partly because of their high star formation rates and, in rarer cases, because 
the submm emission is from an AGN. 

Key words: galaxies: active - galaxies: starburst - galaxies: high-redshift - X-rays: 
galaxies 



1 INTRODUCTION 

There has been vigorous debate about whether the 
most luminous galaxies, most of which have been 
disco vered through their f a r-infrared (FIR) emi ssion 
(e.g. ISanders fc Mirabell Il996l : E owan- Robinsonl [ioooh art- 
violent starbursts, or obscured active galactic nuclei (AGN), 
with the consensu s appearing to be that both processe s 
are important (e.g. IVeilleux et al.|[l995l : iFarrah et al.ll2002h . 
This fits in with the idea of a 'starburst- AGN connection' 
where accretion and star formation, both plausibly triggered 
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by mergers, occur coevally and play a complementary role in 
the formation of galaxies (e.g. ISpringel. Yoshida. fc White! 
l200ll : lHopkins et al.ll2005l ). 

Highly luminous galaxies selected by their submm- 
emission have attracted particular attention. This is in large 
part due to the great observational progress afforded b y 
the SCUBA instrument on the JCMT (|Holland et allll999h . 
which has vastly increased the number of known SMGs. 
While it is clear that su bmm-selected sou rces represent 
a diverse population (e.g. Ilvison et alj l200oh . their nature 
has been the subject of intense study. These efforts stem 
from two key observations: firstly that if SMGs are domi- 
nated by star formation they are among the most power- 
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Table 1. X-ray detected submm sources. Col.(l): Submm source name (|Pope et alj|2005l : IWall. Pope, fc Scotj 120081 ) . Col.(2): Redshift. 
The bold faced redshifts are spectroscopic, otherwise they are photometric. Redshifts in italic are from IRS spectroscopy. Redshifts in 
parentheses are e stimated from IR AC colours, for those sources without optical counterparts. Optical photometric redshifts have ODDS 
parameter > 0.9 llPope et al.ll2005r i, IRAC photometric redshifts have an accuracy of a[Az/(l + z)\ ~ 0.1). Col. (3): 1775 magnitude from 
HST/ ACS imaging. Col. (4): Flux-deboosted 850 fim flux density. Col. (5): IR (8-1000/tm) lum inosity. Ljjj is calculated from the rest frame 
^160 by scaling by the typical Li60-to-LiR conversion from Table A2 of lPope et all [|2006h . Col. (6): Radio luminosity. Col. (7): Chandra 
name Col. (8): Chandra RA. Col. (9): Chandra Dec. Col. (10): Positional offset between IRAC a nd Chandra position (arcs ec) , after correctin g 
for a small (0.6 arcs ec) overall of fset between the two coordinate frames. X- ray d ata are from lNandra et al.l ([2005) and lLaird et all | ||2005| ), 
all other da t a from iPope et al] ||2005| . |2006| , l2008t) . IWall, Pope, fc Scottl [120081 ) and references therein. The redshift for GN10 is from 
iDaddi et ail j2009h . 
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* Included in the A05b sample. The redshift used for GN06 in A05b was 1.865. 
"Southern X-ray and IR counterpart. 
^Northern X-ray and IR counterpart. 



ful starbur sts - and most luminous sources - in the Uni- 
verse (e.g. ISmail. Ivison. fc Blainll 19971 ; iBarger et al1ll998l ; 
iHughes et al] Il998l ; iLillv et al] Il999l ). Secondly, as these 
galaxies are reasonably numerous, they may contribute a 
significant, or even dominant contribution to the global star 
formation rate (e.g. Blain et al. 19991; [Chapman et alj[2005l ; 
lAretxaga et"ai] 120071 ; IWall. Pope, fe Scottll200Sl) . Addition- 
ally, whether or not SMGs are dominated by star formation 
or AGN processes, it seems clear that this sub-mm bright 
phase is a crucial one in the evolution and formation of the 
most massive galaxies. 



Pivotal developments in our understanding of the 
submm population have recently been achieved thanks to 
deep optic al and mid-infrar e d spe c trosco pic, and X-ray ob- 
servations. IChapman et al] ([2003, 2005;) have identified a 
large number of SMGs associated with optically faint radio 
sources in several different fields, following up the counter- 
parts with deep Keck spectroscopy. These observations, in 
addition to providing the first significant sample of secure 
redshifts for SMGs, have shown that they occur at a median 
redshift 2 ~ 2. R ecently, a few rare SMGs have been found 
out to z ~ 4 - 5 (|Capak et alj|2008l ; IDaddi et al]|2009D but 
the bulk of this population is still consistent with a median 
redshift of 2-3. This contrasts with some earlier predictions 



that the SMG population might be at extremely high red- 
shifts. 



Using the IChapman et~aTl (|2005h source list of SMGs, 
Alexander at al. (2005a, hereafter A05a; 2005b hereafter 
A05b) have identified X-ray counterparts to a sample 
of 20 SMGs in the Chandra-Deep Field North (CDF- 
N) region, wher e ultra deep X-ray data are available. The 
IChapman et~ai] (|2005T l sample of SMGs that have spec- 
troscopically identified, optically-faint radio counterparts is 
representative of approximately half of all SMGs. Based 
on the analysis of this sample A05a argued that the AGN 
fraction in SMGs is very high (> 38^^ per cent), and 
hence that the intense star formation in SMGs is accom- 
panied by coeval and near-continuous black hole growth. 
A05b presented a detailed spectral analysis of the X-ray de- 
tected submm galaxies, concluding that 75 per cent of radio- 
selected, spectroscopically identified SMGs host AGN, and 
that the majority (~ 80 per cent) of these are highly ob- 
scured, with column densities > 10 23 cm -2 . Despite the ap- 
parently high AGN fraction, mid-IR spectroscopy of SMGs 
suggests that the AGN is generally no t the dominant con- 
tributor to the bolometric luminosity dValiante et aL I l2007l ; 
IPope et al.l |200S| ; iMenendez-Delmestre et al] I2009T )." so the 
conclusion that they are vigorously forming stars still holds. 

Overall, the above results argue that many of the most 
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massive galaxies in the Universe undergo a period of in- 
tense star formation activity and black hole growth at 
moderate-to-high redshifts, which is deeply shrouded in gas 
and dust. This coeval mode of activity could plausibly re- 
sult in the local corre lations seen between black hole masses 
and bulge properties (iMagprrian et~ln" 1ll998l ; lGebhardt~ et al.l 



l200d ; iFerrarese fc Merrittll200(j). A note of caution should 
be sounded, however, as the " Ichapman et~ail (|2003l , 120051 ) 
method of identifying SMGs, while apparently highly suc- 
cessful, could potentially suffer selection biases due to the 
requirement of radio emission, optical faintness and spectro- 
scopic redshifts. 

Pope et al. (2006; hereafter P06), have used a new re- 
duction of the radio data in the CDF-N together with deep 
Spitzer observations from the GOODS project to identify 
secure counterparts to a large number of high significance 
SCUBA 850 /an sources in the GOODS-N region of the 
CDF-N. The GOODS-N SCUBA survey is small and non- 
uniform but to date it remains the submm field with the 
most complete counterpart identification and deepest X-ray 
data. In this paper, we use the P06 sample to re-examine 
the X-ray emission of SMGs and in particular to assess the 
AGN fraction and X-ray spectral properties. 

Throughout, a standard, flat ACDM cosmology with 
Qa = 0.7 and Ho = 72 km s _1 Mpc~* is assumed. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The Spitzer/ submm sample 

The submm galaxies used in the work have been selected 
from the SCUBA Super-map of the GOODS-N region of 
the CDF-N. The SCUBA Super-map is composed of data 
observed in all SCUBA modes (photometry, jiggle mapping 
and scan mapping) and thus provides inhomoge neous cover- 
age and depth across the GOODS-N field (see iBorvs et al.l 
120031 ; iPope et al.ll2005l and P06 for a full description of the 
Super-map and corresponding analysis) . Taking into account 
the varying depth across the field P06 have presented a 
sample of 35 850 fim selected sources in the GOODS-N re- 
gion, and using radio and Spitzer data have identified highly 
reliable counterparts to 21 of these, with likely counter- 
parts to an additional 12. An update and re-analysis of the 
SCUBA Super-map with additional observation in the fields 
resulted in two further 850 /im selected sources with highly 
reliable counterparts , and is described in the appendix of 
IWall. Pope, fc Scottl (|2008l ). 

For this sample of 35 identified submm galaxies we 
therefore have accurate radio or Spitzer positions, which 
we can then use to unambiguously associate with X-ray 
counterparts. P06 have identified the optical counterparts 
to the SMGs where th ey are detected in the HST images 
l|Giavalisco et al.1 120041 ') . They also assembled together all 
available redshift information, including the available spec- 
troscopic redshifts, photometric redshifts for objects de- 
tected optically, and Spitzer-derived photometric redshifts 
for optically blank SMGs. There are three sources with only 
photometric redshifts and while there is some uncertainty 
in these, particularly in the Spitzer-derived redshifts, for 
our purposes small redshift errors of order Az ~few tenths 
will have little effect on our conclusions. It should be noted 



that GN10 has recently b een identified via CO line emis- 
sion as being at z = 4.042 (|Daddi et al.ll2009l ). considerably 
higher than the previous photometric redshift estimate of 
z = 2.3, and we use the updated redshift value here. Nine of 
the SMGs in t his sample have al so been observed with the 
Spitzer IRS bv lPope et al.l (|2008T ) to confirm their redshifts 
and determine the level of AGN contribution. 

The basic properties of the SMGs with X-ray detections 
are given in Table [T] 



2.2 X-ray data and reduction 

In this paper we use the analysis of the 2-Ms Chandra ACIS- 
I observation of th e CDF-N, which i s the deepest X-ray ex- 
posure yet taken d Alexander et all l2003h . The analysis is 
described in detail in lLaird et al.l (2005). 

Briefly, after applying relatively standard screening pro- 
cedures we created images and exposure maps in several 
bands, notably the full band (0.5-7 keV), hard band (2- 
7keV) and soft ban d (0.5-2 keV). The s ource detection al- 
gorithm described in lNandra et al. I (|2005l ) was then applied, 
which yields a similar, but slightly l arger number of sources 
than used bv l Alexander et al,l (|2003T l. The positions of X-ray 
sources in this catalogue were then correlated with those of 
the submm sources from P06. The positional offsets between 
the X-ray sources and the submm sources, after the removal 
of a small overall offset, are given in Table [T] 

Spectra were extracted for all submm/X -ray associ- 
ated s ources using the methodology described in lLaird et al.l 
(2006). CIAO v3.2.1 and CALDB v3.0.1 were used for the 
spectral extraction and XSPEC vll.3.1 software used for 
the spectral analysis. The energy range of all of the analysis 
was restricted to 0.5-7 keV. Because the CDF-N dataset has 
been taken in several different observing configurations and 
positions, spectra must be extracted separately from each 
observation, and we hav e then summed th ese together using 
the standard FTOOLS (|Blackburnl I1995T ) routines. Source 
counts were extracted from a circular region equal to the 
95 per cent encircled energy fraction (EEF) of the Chandra 
point-spread function (PSF) at 2.5 keV. Local background 
regions were manually selected to avoid contamination by 
nearby X-ray sources. Typically, the background was ex- 
tracted from three to four large regions surrounding the 
source. The effective area and response files were then com- 
bined, weighted by exposure, using the FTOOLS routines 
ADDARF and ADDRMF. As discussed bv lLaird et~aH (|2006h 
the first three Chandra observations were taken at a different 
focal plane temperature to the remainder of the observation. 
We therefore restricted our analysis to the 17 observations 
taken at a focal plane temperature of — 120° C, reducing the 
total exposure of the spectra by 161.7 ks. 

Because of the small number of counts in the individ- 
ual spectra the spectral fitting was p erformed us ing a max- 
imum likelihood method (C-statistic; Cash 1979). This has 
the disadvantage that the g oodness-of-fit cannot be easily 
assessed (but see lLucvll2000h . It is possible, however, to use 
the C-statistic to assess whether a particular model com- 
ponent is needed in a particular dataset and, more impor- 
tantly, to assign confidence intervals to the parameters. All 
of the spectra analysed using the C-statistic were grouped 
into fixed width bins of 4 channels (~ 700 eV), which greatly 
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Figure 1. 2-10 keV X-ray flux vs. 1775 magnitude (AB) for the 
X-ray detected submm galaxies (filled squares). Fluxes are calcu- 
lated for the spectral fits using Model A (see section 3.2). Both of 
the optical/IR and X-ray counterparts to the SMG are plotted for 
GN04 and GN39, and are identified with open squares and circles, 
respectively. The filled circle shows the 2— lOkeV stacking result 
for the undetected SMGs. Not all of the SMGs in the stacked sam- 
ple are detected in the optical, therefore the mean 1775 magnitude 
for the stacked sample is shown as a lower limit. Diagonal lines in- 
dicate constant X-ray to optical flux ratios. The 1775 magnitudes 
were converted using mag= 1.405 X 10 — 7 erg s _1 cm~ 2 . 



improves the processing times in XSPEC without significant 
loss of information for these photon-starved spectra. 

For GN04 and GN39, we identified two potential X-ray 
counterparts (see below). Spectra were extracted from each 
of the individual X-ray counterparts that are associated with 
the IRAC components of the sub-mm galaxies. For GN04, 
where the X-ray sources are only separated by 2.5 arcsec 
and therefore are not fully resolved by the Chandra PSF we 
extracted the spectra using a region smaller than the 95 per 
cent to minimise contamination by the neighbouring source. 



3 RESULTS 

3.1 X-ray and submm associations 

The most basic information is the fraction of submm 
sources associated with X-r ay sources. Cross- correlation of 
the IRAC positions of the iPope et all (|2006l ) SMGs with 
the lNandra et all ([2005) X-ray source list gives a total of 16 
matches within 1 arcsec. The probability that any of these 
are a random association is 0.13 per cent. There are 4 sources 
in the SMGs catalogue which have two IRAC/radio sources 
in the SCUBA error circle, both of which may contribute 
to the submm flux. Of these GN04, GN19 and GN39 have 
X-ray counterparts. In the case of GN19, the X-ray counter- 
part is clea r ly ide ntified with the first counterpart listed by 
IPope et al.l l|200fj| ). i.e. the optically brighter source. In the 
cases of GN04 and GN39 there are X-ray counterparts for 
both of the radio sources. We therefore find an X-ray source 
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Figure 2. Histogram of the effective photon-index of the SMGs, 
assuming no intrinsic absorption (see Table O. The distribution 
of photon-indices for this sample (grey shaded region) and the 
17 sources with effective T in Table 1 of A05b (hatched region) 
are shown. For GN04 and GN39 the photon-indices for both com- 
ponents are shown. The double hatched region shows the A05b 
photon-indices for the 8 sources common to both samples (as in- 
dicated in Table [TJ. Plotting the photon-indices for the common 
sources according to their values in Table[2]from this work shows a 
very similar, but not identical, distribution. Including (excluding) 
the 8 common sources, the K-S probability that the two samples 
are drawn from the same underlying galaxy distribution is 0.02 
(0.005). 



detection fraction for submm-selected galaxies of 46 ± 8 per 
cent, where the error bars are calculated for a binomial dis- 
tribution. 

All the detected submm sources are extremely weak 
in the X-ray, with 2-10 keV fluxes ranging from ranging 
from 4 x 10 -17 erg cm -2 s _1 (close to the limit of the 
survey) to 2 x 10 -15 erg cm -2 s _1 , still a fact or ~ 6 be- 
low the break in the X-ray nu mber counts (e.g. iKim et all 
120071 ; iGeorgakakis et all 120081 ) where most of the X-ray 
background is produced. The X-ray-to-optical flux ratios of 
the detected sources covers a wide range, with many show- 
ing properties in the "normal" range (expected for X-ray 
selected AGN), others relatively optically bright and two 
which are undetected even in the deep HST imaging (7ab > 
28). The Fx/F op t diagram (Fig. [1]) is often used as a diag- 
nosti c tool to help in source classification (e.g.lSchmidt et al.l 
1 19981 ; iGiacconi et~ai1 l200ll ; iHornschemeier et al.ll2003l ). In 
this case, given the extreme faintness of the sources in the X- 
ray, their wide redshift distribution, and the likelihood that 
many of the objects are heavily obscured by dust, these clas- 
sifications are likely to have limited meaning. It is, however, 
clearly evident that SMGs are a heterogeneous class in terms 
of their X-ray to optical flux ratio. 



3.2 X-ray Spectral fitting 
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Table 2. X-ray spectral properties of submm galaxies. All errors correspond to the 90 per cent confidence level. Col.(l): Submm source 
name. Col. (2): Total (source+background) counts in the 0.5-8 keV band. Col. (3): Estimated background counts in source extraction region. 
Col. (4): Hardness ratio (H-S)/(H+S) where H and S are the count rates in the HB and SB, corrected to on-axis values. Col. (5): Exposure- 
weighted off-axis angle. Col. (6): Effective photon index for a fit with only Galactic absorption. Errors Col. (7): Inferred intrinsic JVjj from fits 
to Model A (WABS*ZWABS*(POWERLAW+PEXMON)) with fixed T=1.9 and pexmon normalisation fixed to powerlaw normalisation). 
Col. (8): Improvement in the C-statistic for Model A compared to a T=1.9 power-law spectrum with no intrinsic absorption (Model C). 
Col. (9): Observed frame 2-10 keV flux. Col. (10): Log rest-frame 2—10 keV luminosity (observed). Col.(ll):Log rest- frame 2-10 keV luminosity 
(absorption corrected). Col. (12): X-ray derived SFR, using the lRanalli et aLl ||2003T I relation. Col. (13): Source classification based on X-ray 
properties. 'Amb' denotes sources where the X-ray emission could be from an AGN or star formation. 
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a Spectrum of Southern X-ray source. 
^Spectrum of Northern X-ray source. 



To proceed further in our analysis of these objects re- 
quires a consideration of their spectra. The work of A05b 
has indicated that as a class SMGs can be heavily obscured, 
and some may even be Compton thick. If so, even a quantity 
as simple as the luminosity requires a careful consideration 
of the effects of absorption and hence the spectral shape. In 
particular, at the highest column densities, the AGN may 
only revealed in scattered light, with the spectrum domi- 
nated by Compton reflection. Direct fitting of such models 
to the spectrum can thus constrain the contribution even 
of heavily buried AGN which might be rendered otherwise 
invisible. 

In the modelling of the spectra several different mod- 
els are assumed. To characterise the spectra in the simplest 
manner all of the sources were first fitted by a Galactic 
absorption (Nn = 1.5 x 10 20 cm -2 ) plus power-law model 
(WABS*POWERLAW in XSPEC, Table column 5) to de- 
termine the effective F. In order to characterise any intrin- 
sic absorption we then fit the spectra with two different 
models. The first (Model A) comprises of an intrinsically 
absorbed power- law representing emission from an AGN. 
This characterizes transmission dominated spectra. The sec- 
ond model (Model B) is a neutral reflection model with an 
unabsorbed power-law representing the direct component 
and a reflection component with the same photon index as 



the direct component. This model characterizes reflection- 
dominated spectra, which are often exhibited by Compton 
thick sources. 

Given the very low number of counts in all the objects, 
this is extremely challenging, and this forces us to make a 
number of assumptions during the modeling of the spectra. 
In particular, for the majority of these low-count observa- 
tions it is not possible to constrain simultaneously both the 
power-law photon index, F, and the line of sight absorbing 
column, A^h- For all models we therefore fixed the photon in- 
dex in the power-law and reflection components at F = 1.9. 
The neutral reflection model is represented by the pexmon 
model of Nandra et al. (2007), which include s both con- 
tinuum reflection (Magdziarz & Zdziarski 1995) and the as- 
sociated Fe line emission ( George fc Fa bian 199lj). For this 
model we assumed a single power law illuminating spectrum, 
and the abundances were fixed at the solar value. The as- 
sumed inclination was 60° and the solid angle of the reflec- 
tor Q/2n was fixed to 1. The strength of the reflection is 
therefore parameterized by A TC {, the normalization of the 
reflection component, with the ratio of this to the power 
law normalization, A p i representing the reflection fraction 
R (see Nandra et al. 2007 for details). Note that the val- 
ues of the reflection fraction, R, could be overestimated in 
cases where the soft flux includes a significant contribution 
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from star formation, rather than AGN scattered flux. This 
does not affect the derived maximum AGN luminosity, as 
this depends only on the reflection normalisation (as well 
as the assumed solid angle and inclination of the reflector). 
In model A, we include a reflection component with fixed 
R = 1, which should be a reasonable representation for a 
type 1 (unabsorbed) AGN or a Compton-thin type 2. Model 
B has no intrinsic absorption, but A Te t is free, so the re- 
flection component can dominate, appropriate for Compton 
thick, type 2 sources. The XSPEC representation of Model A 
is therefore WABS*ZWABS(POWERLAW + PEXMON). The 
XSPEC representation of Model B is WABS*(POWERLAW 
+ PEXMON). Despite the additional component in model 
A compared to model B, they have the same number of 
free parameters because in model B the Compton reflec- 
tion strength is allowed to vary. To help assess the improve- 
ment in the spectral fitting due to the intrinsic absorption 
(zwabs) or neutral reflection (pexmon) components we also 
fitted a simple Galactic absorption plus fixed V = 1.9 model, 
WABS*POWERLAW in XSPEC (Model C). For each of the 
three models described above we determine the best fitting 
parameters to the data, as well at the C-statistic value for 
the fit. 

The results of the fits for effective photon index and 
absorbed power-law (Model A) are shown in Table [2] in- 
cluding the C-stat comparison between Model A and Model 
C. The results of the fits to the reflection model (Model B) 
are given in Table [3] A comparison of the C-statistic for the 
three models is given in Table [4] along with the designated 
best fitting model. We adopt a threshold value of AC = 5 
compared to model C. If the improvement in likelihood ex- 
ceeds this value we adopted the more complex model (A or 
B) and choose between A and B based on which delivers 
the smaller C-statistic. The threshold C-value was deter- 
mined using Monte Carlo simulations of power law spectra 
with typical signal-to-noise ratio and background. We per- 
formed 10,000 such simulations, fitting them with an ab- 
sorbed power law (model A) and a reflection model (model 

B) , and comparing the fits to model C . We find that in 1 
per cent of the simulations AC exceeded 5.0 (2.4) for model 
A (model B). The presence of absorption or reflection can 
therefore be inferred at a minimum of 99 per cent confidence 
if AC exceeds this value. 

The distribution of effective photon index for the SMGs 
is shown in Fig. [2] and compared to that from A05b. It can 
be seen clearly that the SMGs in our sample have, on aver- 
age, larger values of F (indicating unobscured X-ray spectra) 
than those from A05. This is supported by results from the 
fits to Model A. As can be seen from Table [2] four of the 
X-ray detected submm galaxies have spectra consistent with 
large intrinsic absorption (Nu > 10 23 cm" 2 : GN04, GN19, 
GN22 and GN39). A further two sources (GN25 and GN40) 
show signs of more moderate absorption in their spectra. 
The majority of the X-ray detected submm galaxies (10/16) 
are therefore consistent with no intrinsic absorption. The re- 
sults of the Model A spectral fitting are consistent with the 
spectral fitting results from A05b for the sources common 
to both samples. 

All 6 sources that show signs of absorption, and only 
those 6, also show large improvements (compared to model 

C) when fitted with model B, the neutral reflection model. 
Results from fits to this model are given in Table [3] 



The reflection fractions range from ~ 5 (GN40), to be- 
ing completely reflection dominated (GN22, GN39-Southern 
source). Comparing models A and B, (Table 2J, we find 
that in GN04 (Northern source), GN19 and GN39 (North- 
ern source) the reflection dominated model gives a lower 
value of the C-statistic, with the other objects showing a 
better fit to the absorbed power law. Typically the differ- 
ence between the two models is not large, however. Excep- 
tions are GN04 (Southern source), GN25, GN39 (Southern 
source) and GN40 where the transmission model is bet- 
ter (all with AC > 5 compared to the other model). In 
GN04 and GN39, which each have two X-ray counterparts, 
each counterpart shows significant obscuration. However in 
each case the Southern source is best fit by the transmission 
model while the Northern source is marginally better fit by 
the reflection model. 

Given the limited photon statistics in our spectra, it 
may not always be possible to identify Compton thick AGN 
purely on the basis of their X-ray spectral shape. This is 
especially so given that the X-ray emission can be heav- 
ily suppressed by large columns of gas. Thus, even relatively 
weak X-ray emitters which do not show evidence for absorp- 
tion or reflection in the spectral fits could still harbour an 
AGN if it is very heavily obscured. It is possible to estimate 
the maximum contribution of a Compton-thick AGN to the 
data using model B, under the assumption that some frac- 
tion of the AGN light is viewed via Compton reflection. We 
adopt the simple assumption that the reflection component 
is seen without any obscuration, whence the normalization 
A Te { should represent the strength of the illuminating power 
law, even if that power law is completely obscured from view. 
As before we assume a slab geometry with an inclination of 
60° and a solid angle Q./2-k — 1, which are typical param- 
eters, but strongly caution that this introduces uncertainty 
as the precise obscuring geometry is not known. The spec- 
tral fits in model B yield an estimate of A Te { and the value 
of the upper limit (calculated using AC = 2.7, correspond- 
ing to the 90 per cent confidence limit) in all cases, even if 
the reflection component does not improve the fit. Table [3] 
shows the best estimate of the intrinsic AGN luminosity de- 
rived using the best estimate of A Te { , as long as it exceeds the 
observed luminosity. We also show the estimated maximum 
AGN luminosity based on the upper limit to A le {, i.e. the 
luminosity of the illuminating power law which would pro- 
duce the maximum reflection component tolerated by the 
data. Again this is quoted only if it exceeds the observed 
luminosity, otherwise the latter is adopted. 



3.3 Stacking of undetected SMGs 

Given the very high X-ray detection rate of the SMGs in 
the sample, it seems likely that the undetected sources do 
have significant X-ray emission that is below the detection 
limit of the CDF-N data. The mean properties of the unde- 
tected SMGs can therefore be constrained by stacking, which 
has proved very useful for other high redshift galaxy popu- 
lations (e.g. iBrandt et a l. 2001; Hornsche meier et al.l l200"ll ; 
iNandra et afll2002l : lLaird et al.ll2005l . 120061 ). 

We use a s tackin g pr ocedure identical t o that describe 
in lLaird et~al"1 (|2005l ) and lLaird et all (120061 ). which is sim- 
ilar to that used by INandra et al. 

I l|2002l ). and we include 
only a very brief outline here. Counts at the positions of the 
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Table 3. Results from the spectral fitting using Model B. Er- 
rors correspond to the 90 per cent confidence level. Col. (1): 
Submm source name. Col. (2): Reflection parameter: the ratio of 
the normalisation parameter between the reflected component to 
the direct component. Col. (3): Improvement in the C-statistic for 
Model B compared to a T=1.9 power-law spectrum with no intrin- 
sic absorption (Model C). Col. (4): Directly viewed 2-10 keV lu- 
minosity from best fit values to model. Col. (5): Intrinsic 2-10 keV 
luminosity of AGN from best fit parameter values for Model B. 
Col. (6): Maximum intrinsic 2-10 keV luminosity of AGN from 
Model B spectral fitting. All luminosities are in erg s — 1 . 



ID 


R 


AC-stat 


log Lx 
obs. 


log Lx 
intr. AGN 


log Lx 
max AGN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


GN01 


< 0.8 


0.1 


43.6 




43.6 


GN04 a 


> 170 


113.8 


42.9 


44.1 


44.2 


b 


17 +5 

lt -6 


31.3 


43.1 


44.1 


44.2 


GN06 


< 15.7 


0.8 


42.5 


43.1 


43.6 


GN10 


< 16.4 


0.8 


42.7 


43.4 


43.8 


GN12 


< 3.0 


0.0 


42.4 




42.9 


GN13 


< 7.2 


0.0 


40.9 




41.8 


GN15 


< 2.2 


0.0 


43.0 




43.4 


GN17 


< 6.2 


0.0 


41.9 




42.7 


GN19 
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37.4 


42.9 


44.1 


44.2 
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12.8 


42.6 


43.8 


44.0 
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0.0 


42.9 




43.5 
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' D -18 


77.6 


42.7 


43.8 


43.9 


GN26 


< 5.9 


0.0 


42.1 




42.9 


GN30 


< 10.5 


0.0 


41.8 




42.8 


GNSg 6 


24 +l3 
z ^-13 


11.2 


42.6 


43.6 


43.8 


a 


> 140 


42.8 


42.5 


43.7 


43.9 


GN40 


3 9+2-1 


13.8 


43.6 


44.1 


44.3 



"Spectrum of Southern X-ray source. 
^Spectrum of Northern X-ray source. 



IRAC counterparts to the SMGs are extracted using a cir- 
cular aperture and summed to find the total counts for the 
SMGs in the sample. These are then compared to the back- 
ground counts, which are estimated from a source-masked 
image using randomly-shuffled positions that are 5-10 arcsec 
from the SMGs. The background estimate is repeated 1000 
times. We use an aperture radius of 1.25 arcsec to extract the 
counts from the soft, hard and full bands. This empirically 
determined value was found to give the strongest signal for 
a sample of stacke d Balmer Break galaxies in the same field 
i|Laird et al.ll2005h and the SMGs in this work. 18 SMGs fall 
within 9 arcmin from the aimpoint of the Chandra observa- 
tions and we include these in our stacking sample. 

The results of the stacking are shown in Table [S] The 
sample of SMGs is significantly detected in the soft, hard 
and full bands, however we strongly caution that the hard 
band signal is sensitive to the extraction radius used and 
the detection significance, a, is only greater than 3 for an 
extraction radius of 1.25 arcsec. Fig. [3]shows the counts de- 
tected in each extraction cell in the soft and hard bands 
and demonstrates that the stacked signal is representative 
of the whole sample and is not dominated by a few sources. 
The fluxes are calculated assuming a spectrum with V = 1.9 
and Galactic absorption. The galaxies in the sample have 
redshifts ranging from 0.8 to 4.1 so to calculate the mean 2- 



Table 4. C-statistic values for models A, B and C. Model 
A: wabs*zwabs*(powerlaw+pexmon) with fixed T=1.9 and pex- 
mon normalisation fixed to powerlaw normalisation. Model 
B: wabs*(powerlaw+pexmon) with fixed T = 1.9. Model C: 
wabs*powerlaw with fixed T = 1.9. 



ID 


C-stat 
Model A 


C-stat 
Model B 


C-stat 
Model C 


Best-fit 
Model 


GN01 


22.4 


21.7 


21.8 
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GN04 a 


33.3 


41.4 
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33.2 


32.6 


63.9 


B 


GN06 


40.4 


39.6 


40.4 


C 


GN10 


31.8 


31.7 


32.5 


C 


GN12 


35.1 


35.0 


35.1 


c 


GN13 


30.4 


30.3 


30.3 


c 


GN15 


28.1 


28.2 


28.2 


c 


GN17 


36.8 


36.8 


36.8 
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GN19 


32.9 


30.4 


67.8 


B 


GN22 


26.1 


28.4 


41.2 


A 


GN23 


40.4 


40.3 


40.3 


C 


GN25 


39.4 


44.6 


122.2 


A 


GN26 


19.6 


19.6 


19.6 


C 


GN30 


37.7 


37.6 


37.6 


C 


GN39 6 


43.9 


41.4 


52.6 


B 


a 


30.8 


39.0 


81.8 


A 


GN40 


32.1 


42.6 


56.4 


A 



"Spectrum of Southern X-ray source. 
''Spectrum of Northern X-ray source. 



10 keV X-ray luminosity of the sample we use a weighted av- 
erage method which assumes that each object in the stack is 
as intrinsically luminous. For every source in the sample we 
calculate the approximate fraction that it contributes to the 
flux by weighting by the ratio of the number of counts in the 
extraction cell for that object to the total number of counts 
detected for the sample. For each object we are then able to 
calculate an approximate luminosity. These are then aver- 
aged to get the mean Lx (2-10 keV) for the sample. In each 
band we convert to rest frame 2-10 keV assuming V — 1.9. 
The X-ray luminosity can be used to provide an estimate 
of the mean star formation rate (SFR) of the sample, un- 
der the assumption that the X-ray emission from the SMGs 
is d ominated by star formation processes and not AGNs 
(e.g. Grimm 1 Gilfanov. fc Sunvaevll2003l ; iRanalli et al.ll2003l ; 



I Persic et al. 2004 ). However it should be noted that the X- 
ray SFR calibration is subject to considerable uncertainty, 
and in galaxies with very high SFRs the X-ray emission may 
be significantly c ontaminated by an obscured AGN (e.g. 
Persic et all 12004) Here we calculate the SFRs using the 



Ranalli et al.1 (|2003) relation; using the IPersic et all (|2004l ) 



relation, for example, would result in SFRs a factor of 5 
lower. 

Using the full band signal, we estimate an average star 
formation rate of 136 ± 24 Mq yr _1 . Interestingly, the soft 
band stack implies a SFR ~ 1.5 times less than the full 
band, at the 1.5-2<r level. This, and the detection of the 
stack in the hard band at all (which at the mean redshift of 
2.1 corresponds to the 6-20 keV rest frame), implies there 
may be significant contamination of the stacked signal by 
AGN. Alternatively, in these very gas and dust- rich envi- 
ronments, star forming X-rays could suffer significant ob- 
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Figure 3. Counts distribution for the undetected SMGs in (a) 
the 0.5-2 keV and (b) 2-7 keV bands. The vertical dashed line 
denotes the mean background counts per extraction cell. 



scuration. However, as noted above, the significance of the 
stacked hard band signal is highly sensitive to the extraction 
radius used and is only > 3<r for the radius used here (1.25 
arcsec). Therefore we hesitate to draw any robust conclu- 
sions from this result. 



4 DISCUSSION 

We have presented a detailed analysis of the X-ray properties 
of submm-selected galaxies in the GOODS-N region. Using 
the 2-Ms Chandra exposure, and accurate positions for the 
SMGs from deep radio and Spitzer imaging, we have been 
able to characterise the X-ray emission of these sources and 
assess their AGN content. For directly detected sources, we 
have used the X-ray spectra to determine accurate luminosi- 
ties and absorption properties. For the undetected sources, 
we have determined the average properties by stacking. 



4.1 Origin of the X-ray emission and AGN 
content of SMGs 

We find that the direct X-ray detection rate of submm- 
selected galaxies in the 2-Ms Chandra data is high, at 46 ± 8 



per cent. The detected sources cover a wide range of red- 
shifts, luminosities and X-ray-to-optical flux ratios. The X- 
ray detected SMGs are thus het erogeneous, much like the 
general SCUBA population (e.g. Ilvison et al.l l200ol ) . Com- 
paring their X-ray luminosities with the FIR and radio lu- 
minosities (Fig. Q, we find that around half of the X-ray 
detected SMGs are consistent with the X-ray emission aris- 
ing purely from stellar processes. Using X-ray luminosity 
as a proxy for SFR (Table , the X-ray derived SFRs for 
these galaxies reach at least 500 Mq yr" 1 , and possibly 
>1000 Mq yr _1 in three extreme cases. GN10 (the highest 
redshift source in the sample) is the most extreme example, 
with an X-ray luminosity of almost 10 43 erg s _1 and star 
formation rate close of ~ 1700 Mq yr . While extreme, X- 
ray derived SFRs are subject to considerable uncertainties, 
espec ially at large SFRs (e.g. iPersic et al.ll20oj ; iTeng et al] 
l2005h . and this is fully consistent with the large SFRs cal- 
culated from the FIR an d polycyclic arom atic hydrocarbon 
(PAH) luminosities (P06. |Pope et al.ll2008h . 

Seven of the SMGs show X-ray luminosities well in ex- 
cess of that expected from star formation, however, and 
these objects almost certainly host AGN. The number of 
unambiguous AGN in our sample is 20 ± 7 per cent, rising 
to 29 ± 8 per cent if we include the three ambiguous cases 
which apparently have very high star formation rates. The 
AGN content of these submm-selected sources is certainly 
very high, but the data do not indicate that strong accretion 
activity is universal in SMGs. The very high X-ray detection 
rate of SMGs is likely to be partially due to their very high 
star formation rates. Furthermore, as we discuss below, at 
least some SMGs are likely to have their FIR flux powered 
by an AGN. It is therefore unclear at this point whether 
the strong overlap between the SMG and X-ray populations 
points towards a connection between intense star formation 
and accretion activity in these sources. 

Our results can be compared with those of A05a, who 
found an even higher AGN fraction than we do, conclud- 
ing that black holes in SMGs are growing near continuously 
during periods of intense star formation. The differing con- 
clusions may seem surprising, as we have studied SMGs in 
the same region (GOODS/CDF-N) as analyzed by A05a, 
and are using essentially the same X-ray data. Indeed 6 (8) 
of the 16 X-ray detected SMGs studied here were also in the 
A05a (A05b) samples (Table [T}. For the common sources, 
our classifcation of the origin of the X-ray emission agrees 
with A05a and A05b, and thus the main differences can be 
attributed to the selection and analysis of the submm sam- 
ples. 

The primary difference is that we are u sing purely 
subm m-selected sources (predominantly from IPope et al.l 
l2005h . whereas A05a use a mixture of SMGs and optically 
faint radio sources with followup SCUBA photometry. The 
AGN fraction of SMGs found in this work (20 - 29 ± 7 
per cent, depending on the origin of the X-ray emission 
in the three ambiguous sources in Table [2} is consistent 
with the lower limit of the AGN fraction quoted by A05a 
(> 38± 

io P er cent). A05a assigned the AGN fraction as a 
lower limit because they were only able to assess the X-ray 
properties of SMGs with radio counterparts and spectro- 
scopic redshifts. The SMGs which were radio-undetected in 
the A05a sample might have been X-ray sources (and hence 
AGN), but without accurate positions they were unable to 
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Table 5. Stacking results of undetected SMGs. All errors are lcr . Col.(l): Observed-frame energy band. Col. (2): Mean redshift. Col. (3): 
Detection significance. Col. (4): X-ray flux per galaxy. 0.5-2keV, 2— lOkeV and 0.5— lOkeV fluxes are given for soft, hard and full bands, 
respectively. Col. (5): 2—10 ke V X-ray luminosity p er galaxy. Luminosities corrected to rest frame 2-10 keV using T = 1.9. Col. (6): X-ray 
derived mean SFR, using the iRanalli et al.l l l2003h relation. Errors are statistical only. 



Band 




17 


Fx 


Lx 


SFR X 








(10~ 17 erg cm" 2 s" 1 ) 


(10 41 erg s" 1 ) 


(Moyr- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Soft 


2.1 


7.84 


1.06 ± 0.20 


4.32 ±0.83 


86 ± 17 


Hard 


2.1 


3.15 


3.18 ± 1.21 


8.22 ±3.13 


164 ± 62 


Full 


2.1 


7.10 


4.31 ± 0.78 


6.81 ± 1.23 


136 ± 24 



tell. As noted by A05a themselves, the radio selection can 
cause a potential bias in favour of selecting AGN. Requir- 
ing a spectroscopic identification can cause a further bias, 
due to the availability of UV/optical lines in AGN spectra. 
Using a deeper radio map, and Spitzer data, we have good 
identifications and accurate positions for the vast majority 
of our SMGs, and we find that the high AGN fraction among 
spectroscopically identified, radio-detected submm galaxies 
fails to extend to the general SMG population. An addi- 
tional difference worthy of note he re is the submm signifi- 
cance threshold. iPope et all {2005) applied a strict limit of 
3.5a to be considered, whereas many of the A05a SMGs have 
lower significance. 

Our excellent positions have enabled us, for the first 
time, to perform a stacking analysis of the X-ray undetected 
submm population. We find a very clear soft band detection, 
and tentative evidence for hard band emission. The mean 
X-ray luminosity suggests an average star formation rate of 
~ 150 Mq yr" 1 , with the full band delivering a higher star 
formation rate than the soft X-ray. Including the 9 SMGs 
from Table [5] that are not classified as AGN dominated then 
the average X-ray derived SFR for the SMGs is ~ 300 M 
yr _1 . This is considerably less than the mean SFR of ~ 
1200 M Q yr" 1 estimated from the FIR luminosity (P06). 
With the strong caveats that calibration of X-ray derived 
SFRs is uncertain (especially at high SFRs) , that the stacked 
signal may be due to emission from obscured AGN, and that 
the SFRs from the FIR may be overestimated, there is thus 
some suggestion that the SMGs may host heavily obscured 
starbursts, with column densities exceeding 10 22 cm -2 . 



4.2 Evidence for obscuration 

A05b presented an analysis of the X-ray spectra of the 
A05a sample, concluding that the majority of X-ray emit- 
ting SMGs are heavily obscured, supporting the interpre- 
tation of the X-ray emission being from an AGN. For our 
(small) sample of 7 SMGs that we classify as AGN we con- 
firm this, and X-ray absorption appears to be an effective 
tool in discriminating between AGN and star forming SMGs. 
All but one (GN01) of the objects we identify as an AGN 
shows evidence for obscuration. The column densities in the 
others are of order ~ 10 23 cm -2 . Three sources are candi- 
date Compton thick objects. GN04 (Northern counterpart), 
GN19 and GN39 (Northern counterpart) have large column 
densities (> 10 23 cm -2 ) and are well fitted by the Compton 
reflection dominated model but the transmission model pro- 
vides almost as good a fit. On the other hand, overall the 



X-ray detected sources in our sample show far less evidence 
for widespread X-ray absorption than was found by A05b. 
From the distributions of effective photon index (Fig. [2J it 
is clear that SMGs in our sample have, on average, larger 
values of F than those in A05b. The K-S probability that the 
two samples are drawn from the same underlying galaxy dis- 
tribution is 0.02. Excluding the 8 sources common to both 
samples, which are selected from the same galaxy popula- 
tion (see Table [T] for details) , the probability that the two 
samples are drawn from the same distribution is 0.005. We 
conclude that the majority (~ 60 ± 12 per cent) of X-ray 
emitting SMGs are unobscured, consistent with the inter- 
pretation that the X-ray emission is from stellar processes. 



4.3 Can AGN activity power the submm 
emission? 

Our X-ray spectral analysis can be used to estimate the 
possible contribution of an AGN to the submm emission. 
For unobscured or Compton thin AGN this is relatively 
straightforward, albeit subject to substantial uncertainty. 
In such cases the absorption-corrected X-ray luminosity in, 
say, the 2-10 keV rest frame band, can be converted using a 
bolometric correction and compared to the FIR luminosity. 
The bolometric correction from 2-10 k eV has been discussed 
exten s ively in the lite r ature fe.g. [Padovani fc Rafanellil 
1 19881 ; lElvis et al l 1 1994 iBarger et all I2005T I. There is a 
very wide range of bolometric correct ions, and suggestions 
that it depends on luminosity (e.g. iMarconi et all 120041 ; 
iHopkins, Richards, fc Hernquisq 12007) and/or Eddington 
ratio (jvasudevan fc Fabianll2007l ). The 'canonical' value of 
Elvis et al (Il994) is W2-10 = 40, but the range given by 



Vasudevan fc Fabianl (|2009h covers approximat ely one or- 



der of magnitude «2-io = 15-150. Adopting the Elvi s et ail 
l| 19941 ) correction, we find that the bolometric luminosity is 
dominated (> 50 per cent contribution) by an AGN in just 
one of the SMGs, that being GN04. If the upper end of the 
possible K2-10 is used then the bolometric luminosities of 
GN19 and GN39 might also be AGN dominated (Fig.©. Ap- 
parently, AGN make little contribution to the total submm 
power in the overall population. 

On the other hand, one could hypothesise - as has been 
argued by e.g. A05b - that many submm sources might har- 
bour very heavily obscured, Compton thick AGN. Using our 
reflection model fits, it has been possible to estimate the 
AGN contribution to the luminosity even in sources where 
we find no evidence for AGN activity. The key point is that 
a luminous AGN should be revealed via a Compton scat- 
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Figure 4. Rest-frame 2-10 keV luminosity for the X -ray detected SMG s against (left) FIR bolometric luminosity and (right) rest-frame 
1.4 GHz radio luminosity. The former are taken from lPope et all d2006l) and the latter have been extrapolated from the observed 1.4 GHz 
fluxes assuming a power law with a = —0.75. The filled circles show the SMGs where the X-ray emission is consistent with that from 
purely stellar processes. The SMGs where the origin of the X-ray emission is either an AGN or ambiguous are shown as filled squares. 
The X-ray luminosities are those calculated using Model A and have been corrected for intrinsic absorption where relevant (Tabic [2}. 
The filled triangles show the stacking result for the undetected SMGs. The 2-10 keV luminosities have been calculated from the soft band 
stacking signal, which was the most significant. The right hand plot shows the mean X-ray luminosity for the 13 SMGs in t he stack with 
radio detections. The solid and dotted lines in both cases show the local relations for purely star forming galaxies from iRanalli et al.l 
j2003h . The dashed line in t he left plot shows the average X-ray to FIR ratio found for a sample of ultra-luminous infrared galaxies that 
are not dominated by AGN llTen g et al. 2005T). For GN04 and GN39 the total FIR, radio and X-ray emission is plotted for these sources, 
as both components are expected to be contributing to the sub-mm flux. The X-ray luminosities of most of the SMGs are consistent 
with this relation, but at high X-ray luminosities (> 10 43 erg s — 1 ) there is a clear departure whereby the objects are dominated by an 
AGN in the X-ray. 



tered X-ray continuum, even when the direct continuum is 
obscured from view. Such a Compton reflection component, 
with accompanying intense iron Ka line emission, is a defin- 
ing charac teristic of Compton thick AGN in the local Uni- 
verse (e.g. lRevnolds et alJll99l : iMatt et alJl200Ch . The lack 
of such a component in the majority of spectra in our sample 
limits the possible contribution of an AGN to the submm 
typica lly to between 1-10 per cent, assuming the 1] vis et al.l 
(1 19941 ) bolometric correction. The exceptions are the three 
objects noted above. Taking the estimates of the maximum 
possible contribution allowed by the spectral fits, and the 
most extreme bolometric correction, we find that GN25 and 
GN40 might also have as much as 50 per cent of their submm 
emission powered by the AGN (Fig. 0. For Compton re- 
flection dominated objects, such estimates are very uncer- 
tain, as the geometry of the reflecting medium (including 
the degree of self-obscuration) is not well known, and this 
will strongly affect the strength of the reflection signatures. 
Overall, however, we can conclude reasonably robustly that 
around 85 per cent of SMGs have their FIR emission domi- 
nated by star formation, whether or not an AGN is present. 

An alternative way of diagnosing the presence of an 
AGN and its contribution to the FIR emission is via m id-IR 
spectroscopy. Valiante et alj 120071 ) , iPope et all l|200Sl ) and 
IMenendez-Delmestre et al] f2009|) have analysed IR spec- 
troscopy of more than 40 SMGs, a few of which overlap with 
the sample analysed here. These works showed that the bolo- 
metric luminosity of most SMGs is starburst dominated and 
that mid-IR dominant AGN in SMGs are rarely found: only 
15-17 p er cent of SMGs are dominated by an AGN in th e 
mid-IR <|Pope et al.ll2008l ; IMenendez-Delmestre et al.ll2009l ). 
For the sources common to this sample, the mid-IR spec- 



tral diagnostics generally agree with the X-ray diagnostics. 
The exceptions being GN19 and GN39 which are both clear 
AGN from the X-ray but are sta rburst dominated based on 
their mid-IR spectral properties (|Pope et alj|200sl ). 



5 SUMMARY 

We have presented an estimate of the AGN fraction and 
X-ray spectral properties of submm-selected galaxies in the 
GOODS-N field, based on the deep (2-Ms) Chandra obser- 
vation. While SMGs show a high X-ray detection fraction 
(45 ± 8 per cent), and also AGN fraction (20 - 29 ± 7 per 
cent), the latter is more modest compared to some previous 
estimates. Most X-ray detected SMGs have low X-ray lu- 
minosity and are unobscured, so are probably dominated by 
stellar emission in the X-ray. On the other hand, the (minor- 
ity) AGN SMGs are obscured with moderate column densi- 
ties, and we find three candidate Compton thick objects. We 
have estimated the AGN contribution to the bolometric lu- 
minosity and find only three objects in our sample in which 
accretion is the dominant power source. The X-ray spectra 
therefore indicate that star formation dominates the FIR 
emission in at least ~85 per cent of the submm population, 
in agreement with mid-IR spectral diagnostics. Overall, we 
conclude that while intense, coeval star formation and in- 
tense accretion is undoubtedly present in some SMGs, this 
behaviour is far from universal, with only around 10-15 per 
cent of the population showing clear evidence for this mode 
of activity. 
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Figure 5. Rest- frame 2-10 keV luminosity for the X-ray detected 
SMGs against FIR bolometric luminosity. The FIR bolometric lu- 
minosities and the symbols are the same as in Fig. [4] The lower 
values of the X-ray luminosities are our best estimate of the in- 
trinsic luminosities for each source, calculated for the best-fitting 
model (Model A, B or C), as given in Table [4] The upper values 
for the X-ray luminosities are the maximum intrinsic values al- 
lowed by the reflection model (Model B, Table |3J and represent 
the maximum possible AGN contribution. The diagonal lines in- 
dicate an AGN-to-FIR bolometric luminosity ratio of 0.5 for 2-10 
keV bolometric corrections of K2— 10 = 15, 40 and 150. For sources 
lying in the shaded region above the lines the total bolometric lu- 
minosity is AGN dominated (> 50 per cent contribution). The 
five AGN where the bolometric luminosity could be AGN domi- 
nated arc identified. 
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